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Abstract. We consider a family of hard core objects moving as independent Brownian 
motions confined to a vessel by reflection. These are subject to gravitational forces modeled 
by drifts. The stationary distribution for the process has many interesting implications, 
including an illustration of the Archimedes' principle. The analysis rests on constructing 
reflecting Brownian motion with drift in a general open connected domain and studying its 
stationary distribution. In dimension two we utilize known results about sphere packing. 



1. Introduction 

We consider a model involving "hard core objects" (typically, spheres) moving as inde- 
pendent Brownian motions, reflecting from each other, and subjected to a constant "force", 
that is, having a constant drift. The objects are confined to a "vessel" by reflection, that 
is, they cannot leave a subset of Euclidean space. Our "toy model" illustrates several well 
known physical phenomena for liquids, under some technical (mathematical) assumptions. 
We prove some theorems for moving objects of any size and shape, but the most interesting 
examples involve a large number of spheres, most of them small. 

The first of the three phenomena that our model generates is tight packing of the objects 
under large pressure, and the formation of the surface of the liquid, that is, a hyperplane 
such that most spheres are below the surface, and there is little room to pack any more 
spheres below the surface. 

The second phenomenon is the "centrifuge" effect. Centrifuges are used to separate materi- 
als consisting of small particles (molecules) with different mass. In this example, we consider 
spheres of the same size but subject to different "forces", that is, drifts. The "heavy" spheres 
tend to be closer to the bottom than light spheres. 

The third phenomenon is Archimedes' principle, which says that an object immersed in 
a fluid is buoyed up by a force equal to the weight of the fluid displaced by the object. We 
will illustrate this principle by a family of two-dimensional discs. One disc is large and it 
is submerged in a "liquid" consisting of a very large number of much smaller discs. The 
large disc either "floats" or "sinks" depending on the ratio of its drift and the drift of small 
discs. This model is limited to the two-dimensional case because the classical sphere packing 
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problem is completely understood only in this case. A similar probabilistic theorem can 
be stated and proved in higher dimensions if the relevant information on sphere packing is 
available. 

Finally, we will give an example involving objects with "inertia" , in which high inertia will 
make the large object (disc) sink more easily, even if the drift of this object would not be 
sufficient to make it sink without inertia. The inertia is modeled by oblique reflection and 
low diffusion constant. Two spheres are said reflect in an oblique way if the amount of push 
(a multiple of the local time) experienced by the spheres during reflection is not identical for 
the two spheres. 

Although Section [2] contains technical material needed for our main results, it may have 
some independent interest. We construct there reflected Brownian motion with drift in an 
arbitrary Euclidean domain and find its stationary distribution. The reason for this great 
level of generality is that we apply these results to the configuration space of hard core 
objects. Even if the objects are spheres, the configuration space does not have to be smooth. 
Some regularity properties of the configuration space for non-overlapping balls were proved 
in Proposition 4.1 in [TP] . 

The present article has its roots in an analogous one-dimensional model studied briefly 
in Section 2 of A construction of an infinite system of reflecting Brownian hard core 
spheres was given in [22]. See the introduction and references in that paper for the history 
and ramifications of the problem. 

We are grateful to Charles Radin for the following remarks and references (but we take 
the responsibility for any inaccuracies). A physical system that could reasonably be called 
a Brownian fluid is a colloid, like milk. Since we are using "reflecting hard spheres", this 
specializes to noncohesive, hard-particle colloids. One of the classic material properties which 
are demonstrated in colloids is the fluid/solid phase transition known by simulation in the 
hard sphere model. Although this was first demonstrated earlier by others, the definitive 
paper seems to be [25J. A short expository introduction to related problems can be found 
online ([23]). A recent preprint concerned with the motion of globules is |15j . 

The rest of the paper is organized as follows. Section [2] is devoted to the analysis of 
reflected Brownian motion with drift in an arbitrary Euclidean domain. Section [3] contains 
some lemmas about the geometry of the configuration space of objects in a vessel. We give 
a sufficient condition for the existence of a stationary distribution for a family of reflecting 
objects in Section HI Finally, Section [5] contains our main results, informally discussed above. 
Our mathematical model is formally introduced at the beginning of Sections [3] and HI 

We are grateful to Amir Dembo, Persi Diaconis, Joel Lebowitz, Charles Radin, Benedetto 
Scoppola and Jason Swanson for very useful advice. 
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2. Stationary distribution for reflected Brownian motion with drift 



Constructing reflected Brownian motion with a constant drift v in a general open con- 
nected set D C ffi" is quite delicate. Even its definition needs a careful formulation. This 
will be done in this section. We will then give a formula for the stationary distribution of 
this process. See [TT1 [121 [19] for some related results. A noteworthy aspect of Theorem 12.31 
below is that we do not require any regularity assumptions on the boundary of D 

When D is a C 3 -smooth domain in W 1 , a (normally) reflecting Brownian motion X with 
constant drift v £ W 1 can be described by the following SDE: 

dX t = dB t + vdt + n(X t )dL t , t > 0, (2.1) 

with the constraint that X t £ D and L is a continuous non- decreasing process that increases 
only when X is on the boundary dD. Here in (12. ip . B is Brownian motion on W 1 and n is 
the unit inward normal vector field of D on dD. When D is C 3 smooth, the strong existence 
and pathwise uniqueness of solution to (12. ip is guaranteed by [20]. When D is a Lipschitz 
domain and v = 0, (12. ip has a unique weak solution by [21 Thm. 1.1 (i)]. Using a Girsanove 
transform, we conclude that weak existence and weak uniqueness hold for (12.11) with non- 
zero constant drift v. When v = 0, this equation has a unique strong solution when the 
domain D is C 7 with 7 > 3/2, by [H Thm. 1.1]. To motivate the definition of reflecting 
Brownian motion with constant drift v in a general, possibly non-smooth, domain, observe 
that when D is C 3 and F x denotes the law of the solution X to (12.11) with X = x, then 
{X,F x ,x £ D} forms a time-homogeneous strong Markov process with state space D. Let 
{Pti t > 0} denote its transition semigroup and p(x) := e v ' x . It is easy to check that for every 
t > 0, Pt is a symmetric contraction operator on L 2 (D,m), where m(dx) := lr,(x)p(x) 2 dx. 
Let (S,^) denote the Dirichlet form of X on L 2 (D,m); that is, 

T = ju £ L 2 (D;m) : limi(u- P t u,u) L 2 (D . m) < 00 j , (2.2) 

£(u,v) = \imj(u - P t u,v) L 2 iD . m) forM,w£j". (2.3) 

Then it is easy to check (see also the proof of Theorem I2.3( i) below) that 

T = {ue L 2 (D;m) : Vn £ L 2 (D;m)} , (2.4) 

S(u,v) — - I Vu(i) • Vv(x)m(dx) for u, v £ J 7 , (2.5) 
2 Jd 

and that (£, JF) is a regular Dirichlet form on L 2 (D; m) in the sense that C C (D) H T is dense 
both in C C (D) with respect to the uniform norm and in T with respect to the Hilbert norm 
a/ Si(u, u) := y/£(u,u) + {u,u)L2(D;m)- Here C C (D) is the space of continuous functions on 
D with compact support. This motivates the following definition. 
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Definition 2.1. Let D C M n be an open connected set. A continuous Markov process 
{Xt, P x , x G -D} taking values in D will be called a reflecting Brownian motion with constant 
drift v if it is symmetric with respect to the measure m(dx) := lD{x)e 2v ' x dx, and ¥ x (Xt G 
dD) = for every x G D and £ > 0, and whose associated Dirichlet form (£, JF) in the sense 
of J22}-(|Z3D is given by «- (l23|) . 

Remark 2.2. (i) Note that reflecting Brownian motion X with constant drift v does not 
have to be a strong Markov process on D. Nevertheless, its associated transition semigroup 
{Pt,t > 0} is well defined by the formula 



J D 

Thus defined {Pt, t > 0} is a strongly continuous semigroup in L 2 (D; m) and so the Dirichlet 
form (SjJ 7 ) is well defined (see [TT]). 

(ii) When v = 0, the Dirichlet form of (E3])-(E1) will be denoted as W 1,2 (Z>)). 

Note that W 1,2 (D) is the classical Sobolev space on D of order (1,2). 

Theorem 2.3. Suppose that v G lR n and D C M n zs open and connected. 

(i) There exists a unique in law reflected Brownian motion X t in D with constant drift v. 

(ii) Suppose that Dk, k > 1, is a sequence of open connected sets with smooth boundaries, 
such that Dk C Dk+i for k > 1 and |Jfe>i = D ■ Let X^ be the reflected Brownian motion 
in with a constant drift v. Assume that Xq — > xq G D in distribution. Then {X^,t > 0} 
converge weakly to {X t ,t > 0}, X = x , in C([0, oo); D), as k — >• oo. 

(iii) The function f(x) = exp(2x-v) ; x G D, is the density of an invariant measure for X . 
If this density is integrable over D then ir(dx) := f(x)dx/ J D f(x)dx is the unique stationary 
distribution for X . 

(iv) Suppose that D and vector vi ^ are fixed, and let v& = 6vi /or 6 > 0. Assume that 
there exists a stationary distribution ir(dx) for reflected Brownian motion in D with drift Vx- 
Then there exists a stationary distribution irb(dx) for reflected Brownian motion in D with 
drift v&, for every b > 1. Suppose that Cq := sup^g^, 2a; • Vx < oo and let D(e) = {x G D : 
2x ■ Vx > c — e} /or £ > 0. Then for every e > we have lim^oo F b (D(e)) = 1. 

Proof, (i) The following facts have been established in [51 E]. Let (^°,JF) be the Dirichlet 
form defined by fl2.4|) - fl2.5p with constant function 1 in place of pi. There is a continuous 
strong Markov process {Y*,¥*,x G D*} on the Martin-Kuramochi compactifictaion D* of 
D, whose associated Dirichlet space is (S°, W l,2 (D)). For 1 < i < m, let /i(a;) := a;,. Then /, 
admits a quasi-continuous extension to D*. Define Y t = (f\(Y t *), • • • , f m (Y t *)), which we call 
symmetric reflecting Brownian motion on D. Process Y admits the following decomposition: 




for every f,g eC c (D). 



Y t = Y + B t + N t 



t > 0, 



M t := exp I / vdB s / |v| 2 cZs ) on each Tt- 
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where B is Brownian motion on M n and N t is an Revalued process locally of zero quadratic 
variation. 

We now construct reflecting Brownian motion X on D with constant drift v through 
Girsanov transform. Let {J~t, t > 0} be the minimal augmented filtration generated by Y*. 
(This notation is similar to the one used for the domain T of the Dirichlet form but there 
will be little opportunity for confusion.) For x G D*, define the measure P x by 

d¥ x ( f* 1 rt 

Since the right hand side forms a martingale under P*, Y* under P x has infinite lifetime. By 
the same (but simpler) argument as that for [7J Lemma 2.4] (with p(x) := e vx there), one 
can show that (Y*,F x ,x G D*) is a symmetric Markov process with respect to the measure 
m(dx) = 1 d{x) p{x) 2 dx . On the other hand, as a special case of [H Theorem 3.1], the 
asymmetric Dirichlet form (£*,JF*) associated with {y*,P.j.} in L 2 (D*\dx) is (£*, W l,2 {D)), 
where 

S*(u,v) = - I Vu(x) ■ Vv{x)dx - [ (v ■ Vu(x))v(x)dx, u,v G W 1 ' 2 {D). 



Denote by {G a ,a > 0} the resolvent of {y*,P x }. The above means that (cf. [211 Theorem 
1.2.13]) 

\u G L 2 (D;dx) : sup(w - /3GpU, u) L 2 {D . dx) < oo \ = W 1,2 (D) 
I P>o J 

and for u,v G W 1 ' 2 (D), 

lim (u — j3G/3U, f)L 2 (D dx) = - / Vn(x) • Vv {x)dx — [ (v • Vu(x))v {x)dx. (2.6) 
P^°° 2 J D J D 

Let (£,T) be the symmetric Dirithlet form of {Y*,¥ x } in L 2 (D*;m). Denote by bW^' 2 {D) 
and bT c the families of bounded functions in W 1,2 (D) and in T with compact support, respec- 
tively. Note that bW^ 2 (D) = bT c , and it is a dense linear subspace in both (W 1 ' 2 (D), (S^) 1/2 ) 
and (F,£{ 12 ). By (1231) . for u G bW^ 2 (D) = bT c , 

lim (u - pGpu, u) L 2 {D . p 2 dx) = lim (u - pGpu, p 2 u) L 2 {D . dx) 

p — >oo p — >oo 

= - I Vm(x) ■ V(p 2 u)(x)dx — j (v • Vu(x))(p 2 u)(x)dx 
2 Jd Jd 

|V 'u(x)\ 2 m(dx) . 

This implies (cf. [17]) that 6JF C c T and for w, v G &.F c , 

£(u, w) = - / Vu(i) • Vv(x)m(<ir) = S(u, v). 
2 
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It follows that T <Z T and £ = £ on T . Conversely, for u G bj- ' c , we have p~ 2 u G bT c . Hence 
hm(u - pGpu, u) L 2 {D . dx) = \im(u - {3G p u, p~ 2 u) L 2 {D . m) 

= £(u,u) 1/2 £(p~\ p- 2 u) 1/2 < oo. 

This implies that u G bW^' 2 (D) C T. In other words, bT c C J? 7 . We conclude that 
(£,J~) = (£1^). This completes the proof that (Y,P X ) is a reflecting Brownian motion with 
constant drift on D. To emphasize that we have constructed reflecting Brownian motion 
with drift as in Definition l2.il we switch to our original notation used in that definition, that 
is, processes Y* and Y under measure P x will be denoted X* and X, respectively. 

Next we establish uniqueness. Suppose that X is another reflecting Brownian motion with 
constant drift v on D. By fT7\ Lemma 1.3.2 and Theorem 1.3.1] the transition semigroup of 
X should be the same as the transition semigroup of X. So as continuous processes, X and 
X share the same law under the initial distribution m. Since the subprocesses of X and X 
killed upon leaving D are Brownian motions in D with constant drift v, it follows that X 
and X have the same distribution for every starting point x G D. 

(ii) Since Dk has smooth boundary, the Martin-Kuramochi compactification of Dk coin- 
cides with the Euclidean closure of D^. So reflecting Brownian motion {Y k ,F*' k } on Dk 
is a strong conservative Markov process with continuous sample paths. Each Y k admits a 
Skorokhod decomposition (cf. [B]) 



Y k = Y k + B k + f n k (Y k )dL k , t > 0, 
Jo 



where B k is Brownian motion on M n , rifc(x) is the unit inward normal vector at x G dD^, and 

L k is the boundary local time for reflecting Brownian motion Y k . As we saw in (i) above, each 

X k can be generated from reflecting Brownian motion Y k on Dk by the Girsanov transform 

rE sk ( f l If* \ 

M k :=exp[ / wdB k -- \v\ 2 ds) on each J*. 



dF* V./n s 2 



x 



Since Yq = Xq is assumed to converge to xq G D in distribution, it is established in j3] that 
(Y k , B k ) converges weakly to (Y, B, P*) in the space C([0, 00), R n x W 1 ) equipped with local 
uniform topology. By the Skorokhod representation theorem (see Theorem 3.1.8 in [13]). we 
can construct (Y k ,B k ) and (Y, B) on the same probability space P) so that (Y k ,B k ) 

converges to (Y,B), P-a.s., on the time interval [0, 00) locally uniformly. Consequently, M k 
converges to M, P-a.s., on the time interval [0, 00) locally uniformly, where 

M t = exp^y v dB s - i J \v\ 2 ds 

Let P be defined by cflP/dP = M t on Tt- Fix T > 0. It suffices to show that X k converges 
weakly to X in the space C([0, T], M. n ). Let $ be a continuous function on C([0, T], M. n ) with 
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< $ < 1. Since $(Y k ) P-a.s., and M| -> M T , P-a.s., by Fatou's lemma, 



Ep[$(F)M T ] < liminfEp[$(F fc )M^] < limsup Ep[$(y fc )M|] (2.7) 



and 



Ep[(l - $)(T)M T ] < liminf E P [(1 - $)(Y*)M£]. (2.8) 

fe— >oo 

Summing (J2ZZD and we obtain Ep[M T ] < limsup^^ Ep[M|]. Note that all M fc 's and 
M are continuous non- negative P-martingales. Hence, Ep[My] = 1 = EpfMp] and, therefore 
the inequalities in (12. 7p and (12.81) are in fact equalities. It follows that 

lim P fc [$(X fc )] = lim Ep[$(F fc )M^] = E P [$(F)M T ] = E P [$(X)]. 

k— >oo fc— >oo 

This proves the weak convergence of X k under P fc to X under P. 

(iii) By definition, m(dx) = l[)(x)e 2v ' x is a symmetrizing measure for reflecting Brownian 
motion X with constant drift v on D. If m(D) < oo, then 7r := m/m(D) is the unique 
stationary distribution of X on D. By Theorem 2(ii) of [16], for every bounded f & 

lim E x [/(X t )] = lim E X [/(X;)] = h(x) for q.e. x E D, 

t— >oo t— >oo 

where /i is a quasi-continuous function with P f /i = h q.e. for every t > 0. Here {P^ t > 0} is 
the transition semigroup of X*. Since P is connected, the reflecting Brownian motion Y* on 
D is irreducible and so is X*. Since m(D) < oo, constant 1 £ J 7 with S(l, 1) = 0. Therefore 
X* is recurrent. It follows that h is constant and equals j D f{x)n{dx). Note that reflecting 
Brownian motion Y* can be defined to start from every point in x G D and has a transition 
density function with respect to the Lebesgue measure in D. As X* can be obtained from 
Y* through Girsanov transform, the same holds for X*. It follows that for every x G D, 

lim M x [f(X t )\ = lim E^/pQ] = / P 1 f(x)n{dx)= [ f{x)v{dx). 

Since is dense in the space of bounded continuous functions on D, the last formula shows 
that 7r is the unique stationary distribution for X* and X. 

(iv) Since 

exp(2c fe) f 

exp(2x ■ Vb)dx < — / exp(2x • vi)dx < oo, 

d exp(2c ) J D 

it follows that n^dx) = exp(2x • \b)dx/ J D exp(2x • vt,)dx is the stationary probability distri- 
bution for reflected Brownian motion in D with drift v^. 

Note that for x G D(e), we have lim^oo exp(— 2(c — e)b) exp(2a; ■ v b ) = oo, so 



lim exp(— 2(c — £)b) / exp(2x • v fe )<ix = oo 

b - > °° JD(e) 
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and, for similar reasons, 



lim exp(— 2(c — e)b) / exp(2x • v b )dx = 0. 

b ^°° JD\D(e) 



It follows that 



Urn n(D \ D ( e)) < lim ^fff' < lim W) ^ " ^ = . 



□ 



3. Configuration space for hard core objects 

In this section, we will start the formal presentation of our model and we will prove two 
lemmas about the configuration space. 

Suppose that d > 2 and D C M. d is open and connected. The set D represents the space 
where hard cord objects may be located. Note that we did not impose any smoothness 
assumptions on dD. 

Consider closed bounded sets Sk C ~R d , k — 1, . . . , N, N > 1. The sets S^'s represent hard 
core objects. We will think about S^s as moving or randomly placed objects, so we will use 
the notation Sk(y) = Sk + y. The diameter of Sk will be denoted by pk- 

The configuration space D C 18L Nd is the set of all x = (a; 1 , . . . ,x N ), x k G M d , such that 
S k (x k ) C D for aU k = 1, . . . , N, and Sj{x j ) n S k (x k ) = 0, for j, k = 1, . . . , N, j ^ k. Note 
that, intuitively speaking, the configuration space should allow for the objects to touch. This 
is not the case of D because we follow the convention in the theory of reflected Brownian 
motion, where the state space is usually represented by an open set, although the process 
visits the boundary of this set from time to time. 

We will prove that the configuration space is connected for two examples of D and S^s. 
We do not aim at a great generality because, first, the problem of characterizing D and SVs 
such that D is connected seems to be very hard, and, second, our main examples in Section 
|5]are concerned with models where connectivity of D is rather easy to see. 

Example 3.1. (i) Suppose that there exists an upper semi-continuous function g : R d_1 — > K. 
such that D = {(21, . . . , Xa) : X\ > g{x2, ■ ■ ■ , Xd)} and all 5Vs are convex., (ii) Suppose that 
D = {(xi, . . . , Xd) : x\ > 0, x\ + • • • + x 2 d < 1} is a one-sided open cylinder, S^s are balls and 
p k < 1 for k = 1, . . . , N. 



Lemma 3.2. If D and Sk's are such as in Example \3.1\ (i) or (ii) then D is pathwise 
connected. 
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Proof. Suppose that x, y £ D and y. We will describe a continuous motion of objects S& 
inside D, such that the initial configuration is represented by x and the terminal configuration 
is y. 

(i) Consider D defined relative to D given in Example 13.11 (i). Let p* = max^ =1 ^ Pk- 
Our argument will involve constants c±, ci > whose values will be chosen later. First, we 

move continuously and simultaneously all objects S& by c\ units in the direction (1, 0, . . . , 0). 
Let z = (z 1 , . . . , z ) denote the new configuration. Next we dilate the configuration by c 2 
units, that is, we fix Si and we move continuously every object S&, k ^ 1, along the line 
segment [z k , z k + (z k — z 1 )^ away from Si, at the speed C2\z 2 — z l \. We move all S&, k ^ 1, 
simultaneously. Note that the objects Sk will not intersect at any time because they are 
convex. Let zi denote the configuration of the objects at the end of the dilation. 

Now we choose c\ and c 2 so large that all objects S& are always inside D, and the distance 
between any two objects is greater than p*, when they are in the configuration zi. 

Next, we start with the configuration y and we use a similar method to move objects 
Sk continuously from configuration y to a configuration u, such that the objects do not 
intersect in the process of moving, they always stay inside D, and the distance between any 
two objects is greater than p* when they are in the configuration u. We make c\ larger, if 
necessary, so that the distance from any object in the configuration Z\ to any object in the 
configuration u is greater than p*. 

At this point, we can move all objects continuously, one by one, from their location in 
configuration Zi to their place in configuration u. We combine motions from x to z 1; then 
to u, and, by reversing an earlier motion, from u to y. 

Consider the set V of all points z 2 representing the locations of objects Sk at all times during 
the motions. The set Y is connected because the motions of the objects were continuous, 
T C D because the objects always stayed in D and never intersected each other, and clearly 
x, y £ T. We have proved that for any x, y £ D, there exists a connected subset of D 
containing both points — this proves that D is pathwise connected. 

(ii) Now consider D defined relative to D given in Example 13.11 (ii). Recall that x = 
far, . . . , x ), and let 2,) represent the center of the fc-th ball. Find a permu- 
tation (tt(1), . . . , tt(AO) of (1, ... , N) such that > x" (2) > > x* {N) . Similarly, let 
y = (y 1 , . . . , y N ) and y h = (y k , . . . , y k ). Let (er(l), . . . , a(N)) be a permutation of (1, ... , N) 
such that yl {1) > y" (2) > ■ > y° x W . Let b = max + 1. Move the vr(l)-th ball 
in a continuous way to a location inside D such that the first coordinate of its center is equal 
to b + 1 and the ball does not intersect the axis of D. Moreover, we move the ball in such a 
way that it does not intersect any other ball or 3D at any time. Next, move the 7r(2)-th ball 
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in a continuous way to a location inside D such that the first coordinate of its center is equal 
to b + 2 and the ball does not intersect the axis of D. We move the ball in such a way that 
it does not intersect any other ball or dD at any time. Continue in this way, until we move 
the 7r(iV)-th ball to a location inside D such that the first coordinate of its center is equal 
to b + iV and the ball does not intersect the axis of D. We move the last ball in such a way 
that it does not intersect any other ball or dD at any time. Such continuous motions are 
possible because we always take the "top" ball from among those remaining in the original 
position, and the diameter of any ball is smaller than the radius of the cylinder D. 

We now move the balls to the configuration y by reversing the steps. First, we move the 
o"(iV)-th ball to the location where its center is y a( - N \ in a continuous way, such that the ball 
does not intersect any other ball or dD at any time. Next, we move the a(N — l)-st ball 
to the location where its center is y°~( N ~ 1 \ in a continuous way, such that the ball does not 
intersect any other ball or dD at any time. We continue in this way until all balls form the 
configuration represented by y. 

Consider the set T of all points z = (z 1 ,. . . ,z N ) representing the centers of all balls at 
all times during the motions. The set T is connected because the motions of the balls were 
continuous, T C D because the balls always stayed in D and never intersected each other, 
and we also have x, y G T. We have proved that for any x, y G D, there exists a connected 
subset of D containing both points — this proves that D is pathwise connected. □ 



Informally speaking, we will assume that all objects Sf. move as independent reflecting 
Brownian motions, with drifts (—0^,0, ... ,0), with Ofe > 0. Formally, the evolving system 
of objects is represented by a stochastic process X t = (X\, . . . , X?) with values in D. In 
other words, the k-th object is represented at time t by Sj.(X]°). We assume that Xf is 
(A r c?)-dimensional reflected Brownian motion in D with drift 



where cif. > for k = 1, . . . , N. 

The m-dimensional volume (Lebesgue measure) of a set A will be denoted \A\ m . 

Lemma 4.1. Assume that D is connected. Let Df, = {x = (x\, . . . ,Xd) G D : X\ = b} and 
a* = min(ai, . . . , ajy). If there is bo 68 such that Dj, = for all b < 6 one? if there is some 
a < a* so that 



4. Existence of stationary distribution 



V = 



((-ai, 0, . . . , 0), . . . , (-ajsr, 0, . . . , 0)) 



(4.1) 



b— >oo 



lim |_D&U-i exp(— 2ab) = 0, 



(4.2) 



then X has a unique stationary distribution. 
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Proof. The uniqueness of the stationary distribution follows from Theorem 12.31 and the fact 
that D is connected. 

In view of Theorem 12.31 (iii), it is enough to show that exp(2x ■ v) is integrable over D. It 
follows from (14. 2 j) that for some c and all b > b , \Db\d-i < cexp(2ab). This implies that 

/ exp(2x • v)dx < / exp(2x • v)dx = 1 / exp(— 2a k x\)dx 
Jn Jd n £J[ Jd 

N poo N 

= ] / \D X1 \d-i exp(— 2cikXi)dxi < J / cexp(2axi) exp(— 2a k xi)dxi < oo. 
k=i J b ° k=i J b ° 

□ 



5. Examples of macroscopic effects 

In this section, we will be concerned with the distribution of the process X t under the 
stationary distribution P, so we will suppress the time variable t, and we will write X for 
X . We will also use the following notation, X = (X 1 , . . . , X N ) and X k = (X k , . . . , X%), for 
k = 1,...,N. 

5.1. Surface of a liquid. 

Theorem 5.1. Suppose that D and Sk's are as in Example \3.1\ (i) or (U), and D satisfies 
the assumptions of Lemma \4-l\ for some a > 0. Fix some a k > for k = 1, . . . , N and let 
c\ = inf xg D Z~2f=i a j x i- Let A > and a& = Xa k for k = 1, . . . , N. Assume that X has 
distribution P. (Note that (X, P) depends on A through drift v in (14.11) . ) 
(i) For any p, 5 > 0, there exists Ao < oo such that for A > A , 

P(a 1 X 1 1 + • • • + a N X? < Cl + 5)>l-p. 

(ii) For any p, 5 > 0, there exists Ao < oo such that for A > Ao, with probability greater than 
1 — p, for every k = 1, . . . , N, for every z E~R d with z\ < —5, we have, 

(s k (x k ) + z)n (d c u\Js 3 (xj)\ ^0. (5.1) 

V j¥* J 

Part (i) of the above theorem says that if the drift of every process X k is sufficiently large 
then the "weighted center of mass" for a typical configuration of S k s is within an arbitrarily 
small number of the infimum of weighted centers of mass over all permissible configurations, 
with arbitrarily large probability. 

Part (ii) of the theorem says that for an arbitrarily small 5 > 0, if the drift of every 
process X k is sufficiently large then with arbitrarily large probability, there is no room in 
the configuration to move any object S k to a new location that would be more than 5 units 
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in the negative rci-direction below the current location of S^. This means, in particular, that 
there are no spherical holes between jS&'s with diameter sup k pk or greater, 5 units below the 
"surface" of SVs, that is, the hyperplane {(xi, . . . , Xd) £ D : X\ — sup fc sup ygSfc ( X fe) y{\. 

Proof, (i) Recall the notation from Theorem 12.31 (iv). We can identify v in (14. ip correspond- 
ing to A = 1 with vi in Theorem 12.31 (iv) . Then 

{x £ D : aix\ H h a N x* < c x + 5} = {x £ D : 2x • vi > c - 26}, (5.2) 

where Cq := sup xeD 2x • Vi. It is now easy to see that part (i) of the theorem follows from 
Theorem 12.31 (iv) . 

(ii) Let ao = in£k=i,...,N a k- Suppose that for some configuration x £ D, there exist k and 
z £ M. d with Z\ < —5, such that (15. ip is not satisfied. Let y represent the configuration 
which is obtained from x by moving Sk from Sk{x k ) to Sk{x k ) + z. Note that y £ D. Since 
y £ D and z\ < —5, we must have a±x\ + • • • + a^Xi > c\ + ao8. We now apply part (i) 
of the theorem to see that the family of configurations x such that (15. 1 p is not satisfied has 
P-probability less than p, if A is sufficiently large. □ 

5.2. Centrifuge effect. In this example, all objects have the same shape but they are 
subject to different forces (drifts ak)- 

Theorem 5.2. Suppose that D and S^'s are as in Example \3.1\ (i) or (ii), Sj = Si for 
j = 2, • ■ ■ , N, and D satisfies the assumptions of Lemma \4-l\ for some a > 0. Assume that 
X has distribution P. Fix some > for k = 1, . . . ,N, let X > and a& = Aa^ for 
k = 1, . . . , N. For any p, 5 > 0, there exists \q < oo such that for A > Ao, with probability 
greater than 1 — p, for every pair j,k £ {1, . . . , N} with ctj > Q-k, we have X[ > X k — 5. 

The theorem says that in the stationary regime, with arbitrarily large probability, if the 
drift is very strong then the identical objects Sk are arranged in an almost monotone order, 
according to the strength of the drift. 

Proof. The idea of the proof is very similar to that of the proof of Theorem I5.1( ii). Recall 
the definition of c\ from Theorem 15.11 Let ao = supja^ — aj : aij > ak} and note that 
ao < 0. Suppose that for some configuration x £ D, there exist j and k such that aj > ak 
and x\ < x\ — 5. Let y £ D represent the configuration which is obtained from x by 
interchanging the positions of Sj{x^) and Sk(x k ). Since y £ D, aj > a^ — a an d x\ < x\ — b~, 
we must have a\x\ + • ■ ■ + a^x^ < c\ + a S. We now apply part (i) of Theorem 15.11 to see 
that the family of configurations x such that aj > ak and x\ < x\ — 8 for some j and k, has 
P-probability less than p, if A is sufficiently large. □ 



ARCHIMEDES' PRINCIPLE FOR BROWNIAN LIQUID 13 

5.3. Archimedes' principle. We will discuss the phenomena of floating and sinking as- 
suming that d = 2 and Sk's are discs. The reason for the limited generality of this example 
is that our argument is based on the classical sphere packing problem. This problem was 
completely solved in two dimensions long time ago (see [TJ]) and it also has been settled 
in three dimensions more recently (see [IB])- The situation is more complicated in higher 
dimensions; see [91 [HI [18] and references therein for details. We will further limit our dis- 
cussion to the cylindrical domain defined in Example I3.1( ii). because this example captures 
the essence of our claims. 

Theorem 5.3. Suppose that d = 2, D is as in Example \3. Si = B(0, 1/2), and Sk = 
B(0,p) for k = 2,...,N, where p < 1/2. Assume that a 2 = a% = ■ ■ ■ = and X has 
distribution P. 

(i) For any p, 5, 7 > 0, there exists p > and N < 00 such that for p G (0,p ) an d 
N > Nq which satisfy the condition p 2 N\/Y2 > 2 — tt/4 there exists ao > such that if 
a 2 > a and a 1 /a 2 := 71 < (vr/ (4vT2) — 7)/p 2 then 

P ( X\ < sup - 1/2 -5] <p. 

\ k=2,...,N J 

(ii) For any p, 5, 7 > 0, there exists p > and N < 00 such that for p G (0, po) an d 
N > Nq which satisfy the condition p 2 N\J~Y2 > 2 — n/A there exists uq > such that if 
a 2 > a and a\ja 2 := 7 2 > (ir/ (4a/T2) + 7) / p 2 then 

P (Xl > 1/2 + 5) < p. 

The theorem is a form of Archimedes' principle. The first part of our result says that if the 
drift of the large ball is smaller than the sum of the drifts of displaced small balls then the 
large ball will "float" , that is, its uppermost point will be at least very close to the "surface" 
of the "liquid" (or above the surface). The second part says that if the drift of the large ball 
is greater than the sum of the drifts of displaced small balls then the large ball will sink to 
the bottom. Both results assume that the system is in the stationary distribution and all 
drifts are large. 

The condition p 2 N\fY2 > 2 — tt/A is needed to make sure that there is an ample supply 
of small discs to make the large disc float. 

Proof, (i) We will use some results about disc packing in the plane from 0, dH dB] • The 
usual honeycomb lattice packing of disjoint discs has density ix j \J~Y2 and this is the highest 
possible disc packing density. 

Consider the unique honeycomb packing S of discs with radii p in the whole plane in which 
some adjacent discs have their centers on the line parallel to the first axis, and one disc is 
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centered at 0. For a set A, we will say that Si, . . . , S& is a honeycomb disc packing in A if it 
contains all discs in S that are contained in A. 

By abuse of notation, we will use | • | to denote the area of a planar set and also the 
cardinality of a finite set. 

Consider arbitrary p,j,5 > 0. We fix some (3 > 2 — 7r/4 and assume that 2 — 7r/4 < 
p 2 N\/Y2 < [3. Note that it is sufficient to prove the theorem for every fixed (3. 

Suppose that for some configuration x G D, we have 

x\ < sup x\ - 1/2 - 5. (5.3) 

k=2,...,N 

Let a\ = 71 and aj~ = 1 for 2 < k < N. Let A = a 2 . Then a±x\ + • • • + a^x^ = 
71X1 + x\ + • • • + . In view of Theorem 15. II (i), it will suffice to show that for small p > 
there exists S\ = 5i(p, N) > such that for some y G D, 

7 ixJ + x\ + ■ • • + xf > llV \ + y\ + • • ■ + y% + 6 X . (5.4) 

If there exist z = (z\, z 2 ) and 2 < k < N such that 21 < — p and 

(s fc (x fe ) + z)n (d c uU =0, (5.5) 

then we choose the smallest k with this property and let 

y = (y\ ...,y N ) = (x\x 2 ,..., x k ~\ x k + z, x k+ \ ...,x N ). 

Then ([51} holds with 5 X = p. 

Suppose from now on that there are no z and k satisfying ( 15. 5ft . Informally speaking, this 
implies that the disc configuration represented by x has a density bounded below by some 
absolute constant c\ > 0. We will now make this assertion precise. Consider a square 

Q = Q(rx,r 2 , p) = {(zi, z 2 ) G M 2 : n - 3p < 21 < n + 3p, r 2 - 3p < z 2 < r 2 + 3p} 

and assume that Q C D and r 2 + 3p < sup fc=2 .jy ^i- We have assumed that there are 
no z and fc satisfying (I5.5p . so B((ri,r 2 ), p) must intersect at least one disc S&(x fc ) with 
2 < k < N . It follows that Sk(x k ) C Q and, therefore, the area of (U2<fc<v ^/c^)) H Q is 
greater than or equal to c\ := 7r/36 times the of Q. 

We have assumed that p 2 Ny/l2 < (3 so sup fc=2 tv x i < f° r some /?! (/?, ci) < 00. 

Let £>i = (0,2/3i) x (-1,1) and L> 2 = £>i \ S^x 1 ). For < 6 < 2/3 l5 define £>* := 
{(xi,x 2 ) G -D 2 : xi < 6}. An upper estimate of the length of dD\ is c 2 := 7r + 4 + 4/3j. Let 
L>3 = D\{p) = {x G D\ : dist(x, dL^) < 2p} and D\ = D b 4 (p) = {x G (£>£) c : dist(x,^) < 
2p}. We have \D\\ < 4pc 2 and \D\\ < 4pc 2 for small p. Hence, the number N\ = Ni(b) of 
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discs of radius p in the honeycomb packing in D\ satisfies 

(\D b 2 \-Apc 2 )/(p 2 Vl2) <N,< (\D b \+Apc 2 )/(p 2 Vl2). 

We cannot pack N± discs in D b 2 \ if \E$\ < \D b 2 \ - 9pc 2 . If p is small then this condition 
follows from b — b\ > 9pc 2 . So in any configuration of N% discs in D 2 , the iVi-th disc from 
the bottom will be at most 9pc 2 units below the position of the TVj-th disc from the bottom 
in the honeycomb packing of D 2 . 

If & 2 = & + pv^3 then the honeycomb packings of D 2 and D 2 2 differ by one row of discs, or 
a part of one row. The centers of the discs in the top row of D 2 2 are p\f?> units above the 
centers of discs in the top row of D 2 . There are no more than 1/p discs in the top row of 
D 2 2 . Consider any N3 such that N\(b) < N3 < N\(b 2 ). For any configuration of N$ discs in 
_D 2 , the iV3-th particle from the bottom will be at most 9pc 2 + pV3 units below the position 
of the iV 3 -th particle from the bottom in the honeycomb packing of D b 2 . 

Let S 2 (z 2 ), . . . , Sn{z n ) be the disc configuration in D 2 obtained by taking N — 1 discs in 
the honeycomb packing of D 2 with the lowest first coordinates. Then 

x\ > z\ - (9c 2 + V$)p for 2 < k < N. (5.6) 

Therefore for any 2 < N 2 < N, 

x\ + --- + x^>zl + --- + z^- (N 2 - l)(9c 2 + V3)p. (5.7) 

Recall that we assume that there do not exist z and k satisfying (15.51) . Suppose that 
sup fc=2 j7V x\ > sup k=2 N z^ + 5/2. Let /Ci = {1 < k < N : x\ > sup k=j ^ N z{ + 5/2}. 
Note that |/Ci| > Ci5/(2np 2 ) for small p because x k, s represent a configuration with density 
bounded below by c±. We have for small p, 

Let = N — |/Ci|. We apply (15.71) to the discs corresponding to 1 < k < N, k ^ /Ci, to 
obtain 

4- ^>-(N 2 -l)(9c 2 + V3)p. 

2<k<N,k^K.i 2<k<N,k£Ki 

Combining both estimates, we see that, for small p, 

x\- J2 4>Ci5 2 /(8np 2 )-(N 2 -l)(9c 2 + V3)p 

2<k<N 2<k<N 

> Cl 5 2 /(87rp 2 ) - {(3/{p 2 VV2) - l)(9c 2 + y/3)p 

> Cl 5 2 /(167ip 2 ). 
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Hence, for small p > 0, (15.41) holds with y 1 = x 1 , y k = z k for 2 < k < N, and 8\ = 
ciS 2 /(167rp 2 ). Note that y G D because S 2 (z 2 ),..., Sn(z n ) is a part of the honeycomb 
packing of D 2 , so these discs are disjoint and they are disjoint with Si (x 1 ). 
Next suppose that 

sup x\ < sup z\ + 5/2. (5.8) 

k=2,...,N k=2,...,N 

Recall 7 from the statement of the theorem. We can assume without loss of generality that 
7 G (0, 1). If p is small then we can find y\ G (x\ + 5(1 — 7)/4, x\ + 5/4) such that the line 
M := {(ui,u 2 ) : U\ = (x\ + y\) /2} is a line of symmetry for the packing S. Let y 1 = (y\,xl). 
Note that Sii^y 1 ) \ S^x 1 ) is "filled" with the discs from the family S^z 2 ), . . . , Sn(z n ) when 
p is small because, in view of (15.31) and (15.81) . 

sup z\ > sup x k - 5/2 > x\ + 1/2 + 5 - 5/2 > y\ + 1/2 + 5/4. 

k=2,...,N k=2,...,N 

Let : R 2 -»• M 2 be the symmetry with respect to M. Let IC 2 = {2 < k < N : S k {x k ) fl 
Si(y 1 ) 7^ 0}. For /c G /C2, let y fc = 7W(x fe ). For all other x fc , let y k = x k . It is easy to see 
that y G D. 

The next part of our argument is best explained using physical intuition. Suppose that 
all discs Sfc are made of material with mass density 1 and they are in gravitational field 
with constant acceleration 1 in the negative direction along the first axis. When we move 
disc S\(x l ) to the new position at Sid/ 1 ) then we do (y\ — x\)tx/4 units of work, which 
is at least (5(1 — j) / 4)tt / 4. We can imagine that the mass in S^x 1 ) fl S^y 1 ) does not 
move, and we only move the mass in Si(x 1 ) \ S^y 1 ) to its symmetric image Si(y 1 ) \ Si(x 1 ) 
under Ai. When p is very small, the discs Sk(x k ), k G JC 2 , have total mass arbitrarily close 
to (?t/-\/T2)|5'i(?/ 1 ) \ Si (a; 1 ) I, uniformly spread over S^y 1 ) \ S^x 1 ). Hence, the amount of 
work needed to move all discs Sk(x k ) to Sk(y k ) for k G JC 2 is negative and smaller than 
— (ix / \/Y2)(5(l — 2 r y)/4)n/4, for small p. In other words, 

iVi - 4)*P 2 < -(tt/v^X^I - 2 7 )/4)tt/4. 

Recall the assumption that 71 < (ir/ (4vT2) — 7)/p 2 . We have y\ < x\ + 5/4 so 
7i(2/i " x\) < 7 i5/4 < ((n/(4Vl2) - 1 )/p 2 )5/4. 
Combining the last two estimates we obtain 

keK.2 

< -(k/VU)(6(1 - 2 7 )/4)(vr/4)/(7rp 2 ) + ((^(4^) - l)/p 2 )5/4 

< -i5/(12p 2 ). 
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Thus, if we take 5 1 = 7<5/(24p 2 ), then 

7ix} + x\ H V xf = 7ixJ + x\ + x{ 

= + Yl x i + S ^ 

fc£/C 2 \<k<N,k£K.<2 

> iiy\ +J2y k i + 6i+ Yl 

fcG/C 2 l<k<N,k(K 2 

Hence, condition (15.41) is satisfied. This completes the proof of part (i) of the theorem. 

(ii) The second part of the theorem can be proved just like the first part. The proof 
is identical up to (15.81) . In the part following (15. 8p . all we have to do is to take y\ £ 
(x\ — 5/4,x\ — 5(1 — 7)/4) instead of y\ £ (x\ + 5(1 — j)/4:,x\ + 5/4), because in this part 
we want to move Si(x 1 ) down, not up. We leave the details to the reader. □ 

5.4. Inert objects. We will model inertia of objects Sk by changing the rules of reflection. 
When two different objects Sj and Sk reflect from each other, they will no longer receive 
the same amount of push to keep them apart. One way to formalize this idea is to say that 
when the process X hits the boundary of D at the time when Sj hits Sk, then X is not 
reflected normally but it is subject to oblique reflection. The drift will not be assumed to 
be related to the value of inertia for Sk- In other words, a^'s may model forces which have 
strength dependent on factors other than the inertial mass. 

We will assume that the standard deviation for oscillations of Sk is inversely proportional 
to the inertia of Sk- The reason for this assumption is purely technical. The assumption 
allows us to transform the problem to the model covered by Theorem 12.31 In general, it 
is not easy to find an explicit formula for the stationary distribution of reflected Brownian 
motion with oblique reflection (even if the process has no drift). 

Next, we formalize the above ideas. Let rrik > be the parameter representing the inertia 
for Sk- Let 

T(x) = (m 1 x l ,...,m N x N ), x £ R Nd , 
D = T(D), X t = T(X) t . 

We assume that X is reflected Brownian motion in D (with the normal reflection), with drift 

v = ((-a 1 m 1 ,0, ... ,0), ... , (-a N m N ,0, . . . ,0)). 

Let P denote the stationary distribution for X and let P' be the corresponding stationary 
distribution for X. Note that under P', the quadratic variation process for X k is t/m 2 k . 
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In the present example, if two objects Sj and reflect from each other then the infini- 
tesimal displacement of Sj is mk/rrij times the infinitesimal displacement of Sk- 
We will illustrate the effect of inertia using the same model as in the Section 15.31 

Theorem 5.4. Suppose that d = 2, D is as in Example \S '. l\( ii) . Si = B(0, 1/2), and Sk = 
B(0,p) for k = 2,...,N, where p < 1/2. Assume that a 2 — °3 = • " " = a N, m i > 1; 
m 2 = • • • = tun — 1; and X has distribution P'. 

(%) For any p, 5, 7 > 0, there exists po > and N < 00 such that for p £ (0,po) an d 
N > N which satisfy the condition p 2 N\/Y2 > 2 — n/A there exists a® > such that if 
a 2 > a and a 1 /a 2 := 71 < (n/ (4vT2) — 7)/ (p 2 mi) then 

P' [Xl < sup X^ - 1/2 -5) <p. 

V k=2,...,N J 

(ii) For any p, 6, 7 > 0, i/iere exists po > and N < 00 snc/i i/iai /or p £ (0,po) 
N > N which satisfy the condition p 2 N\/Y2 > 2 — n/A there exists Oq > such that if 
a 2 > ao and a\ja 2 := 72 > (tt/ (4 vT2) + 7) / (p 2 m 1 ) i/ien 

P' (Xl > 1/2 + 5) <p. 

The theorem says that the higher is inertia mi, the lower is the critical drift a± that makes 
the disc S\ sink. We note parenthetically that behavior of real particulate matter can be 
paradoxical, unlike in our example. Large and heavy particles may move to the top of a 
mixture of small and large particles under some circumstances; see [24] . 

Proof. We can use the same reasoning as in the proof of Theorem 15.31 but with a twist. 
Theorem 15. II (i) must be applied to the process X under P, so we have to analyze miaix\ + 
a 2 x\ + - ■ - + a^Xi rather than aix\+a 2 xf + - ■ - + a]\fXi . Hence, 71 in H5.4[) must have an extra 
factor of ljm\. The constant 71 in the present theorem has that extra factor, as compared 
to the constant 71 in Theorem 15.31 With this change, the proof of Theorem 15.31 applies in 
the present context. □ 
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